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The momentum and energy of phonons in a Bose-Einstein condensate are measured directly from 
a time-of-flight image by computerized tomography. We find that the same atoms that carry the 
momentum of the excitation also carry the excitation energy. The measured energy is in agree- 
ment with the Bogoliubov spectrum. Hydrodynamic simulations are performed which confirm our 
observation. 



In a Bose-Einstein condensate (BEC), the long wave- 
length excitations are phonons, characterized by a linear 
dispersion relation, as well as a wave function consisting 
of atoms with positive and negative momenta. Both the 
dispersion relation and the wave function were measured 
indirectly, by Bragg spectroscopy Q || . 

A two-photon Bragg process coherently imparts an en- 
ergy ftto and momentum hk, determined by the frequency 
difference and the angle between the Bragg beams, to the 
condensate ||. In Bragg spectroscopy, the two-photon 
energy hio is varied, and the response of the BEC to the 
two- photon Bragg transition is observed JI| || [0. Thus, 
the directly measured energy is the photon energy. Cal- 
culating the excitation energy requires knowledge of the 
process by which the photons impart momentum and en- 
ergy to the BEC. Analogously, excitations in superfluid 
4 He were observed indirectly, by measuring the momenta 
and the energies of scattered neutrons k». 




FIG. 1. (a) Absorption picture of a BEC excited with 
fc£ = 0.96 after 38 msec TOF. The left cloud corresponds to 
the condensate and the right cloud to the released phonons. 
(b) The density distribution of (a) reconstructed by comput- 
erized tomography. A clear separation is evident between 
the two clouds. The energy is computed within ellipsoids of 
sizes indicated by the ellipses. The area of each image is 
400/im x 400^m. 

The phonon energy is much larger than that of a free 
particle with the same momentum, and consists mostly 
of interaction energy between the phonon and the con- 
densate. When the trapping potential is turned off, all 
the interaction energy is transformed into kinetic energy 



during a short acceleration period j?J . 

In this letter we present a direct measurement of this 
energy from time-of-flight (TOF) images of the released 
atoms. Previously the ground state interaction energy 
of the condensate was also measured from TOF absorp- 
tion images [Q. Since absorption images provide only 
the density of the cloud integrated along the axis of the 
absorption beam, whereas energy measurements require 
knowledge of the full density distribution, fits to model 
functions of the density were needed 0. Here we use 
computerized tomography of the TOF absorption images, 
in which the cylindrical symmetry of the cloud is used to 
reconstruct the full density profile from a single absorp- 
tion image. 

Several measurements are performed using this tech- 
nique. First, the ground state interaction energy of a 
released, unexcited BEC is measured and compared to 
the value extracted from a fit to the radial size of the 
absorption image. Second, the density distribution of a 
Bragg-excited condensate is reconstructed. A clear sepa- 
ration between a released-phonon cloud and a condensate 
cloud, which is not clear in the original absorption image, 
is visible. The energy and momentum added to the cloud 
by the Bragg pulse are measured from the reconstructed 
image. In this way the excitation energy is found without 
sweeping through frequencies, and by looking directly at 
images of the atoms. Finally, the energy carried by the 
released phonons is measured by observing the distribu- 
tion of the total energy between the two clouds. Thus, 
we extract the excitation energy from a single absorption 
image. Hydrodynamic simulations of the phonon release 
process are performed and found to be consistent with 
our observations. 

In a previous experiment the interaction energy be- 
tween two identical condensates, with a large relative ve- 
locity, was measured by their interference pattern . For 
this case a one-dimensional treatment was sufficient. 

The relation between E^, the energy, and hk, the 
momentum, of an excitation in a homogenous conden- 
sate has the Bogoliubov form Ek = y/ E r (E r + 2gn) || , 
where n is the condensate density, g is the coupling con- 
stant, given by g — 4irh 2 a/m, m is the mass of the atoms, 
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a is the scattering length, and E r = (Hk) /2m is the cor- 
responding free-particle energy. According to the Feyn- 
man relation, the probability of creating an excitation 
with momentum hk is proportional to E r /E k [ ^0| . Exci- 
tations with k < , where £ = y/H 2 /2mgn is the heal- 
ing length, are called phonons whose dispersion relation 
is approximated by the linear relation E k « gn/mk 

0- 

For an inhomogeneous condensate, we use the local 
density approximation (LDA) according to which the sys- 
tem behaves locally as a part of a homogenous conden- 
sate. The excitation energy is then defined by the average 
of the local excitation energies over the entire condensate 
0, 
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where e k (t) is the excitation energy corresponding to 
1 E 

n(r), S k = — J — -j-^n{r)dr is the structure factor of 

the condensate, and N is the number of atoms. 

In our experiment a nearly pure BEC of approximately 
10 5 87 Rb atoms in the \F,m,f) — |2, 2) ground state, is 
formed in a QUIC type magnetic trap Q. The thermal 
component of the cloud is measured to be 5% or less. 
The trap is elongated and cylindrically symmetric, with 
radial and axial trapping frequencies of 220 Hz and 25 Hz, 
respectively. The radial and axial Thomas-Fermi (TF) 
radii of the condensate are thus Rtf = 3 yum and Ztf = 
27 fj,m |l4} |. £ _1 = 4.3 /im~ 1 averaged over the entire 
condensate. 

The two Bragg beams, produced by the same laser, are 
detuned 6.5 GHz from the 5S1/2, F=2 — ► 5P3/2, F'=3 
transition. The beams pass through two acousto-optic 
modulators that control the frequency difference between 
them. Bragg pulses of duration varying between 1 and 3 
msec are applied to the condensate. The angle between 
the two beams is varied to produce excitations with vary- 
ing momenta along the z axis (to ensure cylindrical sym- 
metry). The beams intensities are adjusted such that 
the number of atoms in the excitation is always less than 
25% of the total number of atoms in the condensate. The 
frequency difference between the two beams is varied to 
produce excitations of different strengths. In the range 
of fc's measured, the dominant broadening of the Bragg 
resonance is due to the inhomogeneous density in the 
condensate ]l2| , and is always narrower than the finite- 
time broadening of our pulses. Since the pulse frequency 
spectrum is broad compared to the excitation line width, 
excitations are formed at every point in the condensate 
JT5| . Thus, the excitation energy is given, according to 
Eq.(l), by the average energy added per excitation, i.e. 
the energy added to the cloud by the pulse, divided by 
the number of excitations. 

After the Bragg pulse, the magnetic trapping potential 
is rapidly shut off, and repulsion between atoms causes a 



short acceleration period, during which all the interaction 
energy is transformed into kinetic energy. After 38 msec 
of free expansion the cloud is imaged by an on-resonance 
absorption beam, perpendicular to z. Figure la shows 
an absorption image for excitations with k^ — 0.96. Two 
overlapping clouds are seen. The left cloud corresponds 
to the condensate, and the right cloud to the released 
phonons. The right cloud is clearly larger in the radial 
direction than the left cloud, reflecting the interaction 
energy of the phonons. 

In the TOF image, the atom's position is proportional 
to their momentum. Therefore, if the density distribu- 
tion n(r, z) after free expansion is known, we can directly 
calculate the energy and momentum of the atoms. The 
absorption picture only provides n(r, z) integrated along 
the absorption beam axis f(y, z) = J n(r, z)dx. However, 
since the cloud is cylindrically symmetric we reconstruct 
n(r, z) from a single absorption image using computer- 
ized tomography. According to the Fourier slice theorem, 
n(r, z) is reconstructed by evaluating fl^ j 

n(r ' z ) = 77TT2 / F ( K v> z ) J o ( K y r ) dK v ( 2 ) 
(2tt) J 

where F(k v ,z) is the ID Fourier transform of f(y,z) 
along the y direction, and Jo is the zeroth order Bessel 
function. The center of coordinates (r = Q, z = 0) is de- 
termined by a Gaussian fit to the BEC cloud. Figure lb 
shows the reconstructed n(r, z) of Fig. la. The released- 
phonon cloud is seen to have a shell-like shape. A clear 
separation between the two clouds, which is not visible 
in Fig. la, is now seen. This separation is caused by the 
released repulsion energy between the two clouds. 

From the reconstructed n(r, z), the energy of the atoms 
is calculated as, 
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where t = 38 msec is the time of flight, and Zi n (z) is the 
in-trap z position of an atom with position z after TOF 
(ji n is approximately 0). For the BEC cloud we take 
Zin{z) = {z/Z final) %TF, where Z fmal is the parabolic 
radius in the z direction of an integrated parabola, fit- 
ted to the absorption image . For the released-phonon 
cloud, we have no theory for Zi n (z). Thus, we take 

Zin{z) = 0. 

First, the energy of an unperturbed BEC is mea- 
sured. To reduce background we measure the BEC 
energy within an ellipsoid with radii obtained from a 
parabolic fit to the absorption image. The average en- 
ergy per atom is measured in 12 images of unperturbed 
BEC's to be EbecI^ 1 — 640 Hz, which is approximately 
equal to the value of EbecI^ 1 — 590 Hz, obtained by 
the more conventional technique of fitting an integrated 
parabola to the same absorption images. 

The energy of a Bragg-excited BEC is then measured 
in the combined volume of two ellipsoids, centered around 
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the BEC cloud, and the released-phonon cloud (see Fig. 
lb). The size of the second ellipsoid is enlarged to in- 
clude all the atoms in the released-phonon cloud. A ref- 
erence image with no excitations is subtracted in the area 
of the released-phonon cloud to exclude the energy car- 
ried by the thermal atoms. The measured momentum 
of the cloud, divided by hk, yields the number of exci- 
tations. During the pulse the average energy per atom 
in the condensate cloud Ebec decreases by ~12%. We 
take the time-averaged value Ebec I 'h — 600 Hz to be 
the relevant value. 
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FIG. 2. The filled circles are the measured energy in the 
condensate and released-phonon clouds as a function of the 
number of excitations for fc£ = 0.96. The open circles are 
the energy measured in the released-phonon cloud only. Solid 
lines are the linear fit to the data. 



of the atoms in both clouds as a function of the mea- 
sured number of excitations for k£ = 0.96. The slope of 
these points is E^/h — 1675 ± 290 Hz and corresponds 
to the energy of a single excitation. It is in agreement 
with Ek/h — 1780 Hz, calculated from the Bogoliubov 
dispersion relation in the LDA and also with the mea- 
sured value of @ E k /h = 1750 ± 20 Hz. Figure 3 shows 
the excitation energies measured this way at four differ- 
ent momenta. The measured energies agree well with the 
Bogoliubov excitation spectrum for our measured Ebec 
in the LDA with no fit parameters, indicated by a solid 
line. 
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FIG. 4. The linear density of the condensate cloud 
(solid line) and the released-phonon cloud (dashed line), for 
fcf; = 0.96 after 38 msec of flight, (a) z direction, (b) r direc- 
tion. The plot shown is an average of several images. The 
position of the Bragg recoil from the center of the condensate 
is indicated by the vertical dotted line. The released-phonon 
cloud is clearly larger than the condensate cloud in the ra- 
dial direction. The condensate cloud is not distorted by the 
release process. 
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FIG. 3. Measured excitation energies, for excitations of 
various k extracted from the slope of the energy in both clouds 
versus the number of excitations. The solid line is the Bogoli- 
ubov dispersion relation, calculated in the LDA for the mea- 
sured Ebec with no fit parameters. The dashed line shows 
the free-particle parabolic dispersion. 

The filled circles in Fig. 2 are the measured energy 



In the reconstructed n(r, z) of Fig. lb, the condensate 
cloud and the released-phonon cloud are distinct and sep- 
arate. Figure 4a and 4b show the linear density profile 
in each of the clouds, along the z and the r directions 
respectively, for kt; = 0.96. The released-phonon cloud 
is significantly larger than the condensate cloud in the r 
direction. 

The condensate cloud seems not to be distorted by the 
interaction with the phonons. This may imply that each 
phonon is repelled by a fraction of the condensate that 
is much heavier than the phonon. Since the condensate 
cloud is much heavier than the phonon cloud, this would 
suggest that most of the excitation energy is carried by 
the released-phonon cloud. In order to check this hypoth- 
esis, we measure the energy of each cloud separately. The 
open circles in Fig. 2 are the measured energy in the 
released-phonon cloud only, as a function of the mea- 
sured number of excitations. From a linear fit, the slope 
of these points is 2390 ± 90 Hz. When the average en- 
ergy of an atom in the condensate cloud is subtracted 
from this value [0, we get an energy of 1840 ± 100 
Hz per excitation in the released-phonon cloud, which is 
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consistent with the Bogoliubov value E^/h = 1780 Hz. 
This indicates that all of the excitation energy is indeed 
carried by the atoms in the released-phonon cloud. 

The filled circles in Fig. 5 are the excitation energies, 
measured from the slope of the energy of the released- 
phonon cloud versus the number of excitations. Except 
for the fc£ = 0.32 point, the measured points agree with 
the Bogoliubov spectrum in the LDA, indicated by a solid 
line. 




FIG. 5. Excitation energies measured from the re- 
leased-phonon cloud only. The filled circles are measured from 
the slope of plots like the open circles in Fig. 2. The open 
circles are measured from a single image. The solid line is 
the Bogoliubov dispersion relation, calculated in the LDA for 
the measured Ebec with no fit parameters. The dotted line 
is the average energy of the released-phonon cloud minus the 
condensate average energy, by a hydrodynamic simulation. 

Given that all the excitation energy is carried by the 
released-phonon cloud, this energy can be extracted from 
a single image. This technique requires the collection of 
far less data than does the indirect measurement, where 
an entire spectrum is needed. The open circles in Fig. 
5 are the excitation energies extracted from the image 
with the largest number of excitations, for each k. These 
points agree with the filled circles measured from the 
slope. 

During the release process some of the interaction en- 
ergy is transformed into motion in the z direction, caus- 
ing the condensate to recoil. The zero of momentum 
is therefore not at the center of the condensate cloud, 
introducing a systematic error in our procedure. For 
the higher k values the interaction energy is transformed 
mainly into motion in the radial direction so this system- 
atic error is negligible. For the smaller k values a larger 
fraction of the energy is transformed into motion in the 
z direction. This may explain the low energies measured 
at the fc£ = 0.32 point. 

We perform hydrodynamic simulations of the release 
of the excited cloud |I3]. Excitations at different initial 
positions in the condensate are chosen with a probability 

— T-r-n (r) . The condensate is then released and the ex- 



citations propagate with an initial momentum hk. The 
force at each time step is equal to the gradient of the 
interaction energy in the condensate at the current po- 
sition of the excitation, multiplied by a factor, account- 
ing for the excitation's initial interaction energy. Since 
the excitations are a small perturbation, the condensate 
is assumed to expand as would an unperturbed conden- 
sate [jjj . The calculated average energy of the excitations 
minus the average energy of the BEC, after 38 msec of 
free expansion, is shown by the dotted line in Fig. 5. 
The simulations are seen to be consistent with the obser- 
vation that the released-phonon cloud carries all of the 
excitation energy. 

In conclusion, the interaction energy of phonons in a 
Bose-Einstein condensate is directly measured by com- 
puterized tomography of TOF images. The excitation 
energy is extracted from a single absorption image. This 
measurement of the excitation energy relies solely on im- 
ages of the atoms, and does not require any information 
about the photon energy. Furthermore, no assumptions 
are made with respect to the shape of the cloud, except 
for cylindrical symmetry. The measured energy is found 
to agree with the predicted Bogoliubov spectrum in the 
LDA. After the cloud is released, the excitation energy 
is seen to be carried by the same atoms that carry the 
excitation momentum. Hydrodynamic simulations of the 
release process are performed, and are found to be con- 
sistent with our experimental findings. 

Computerized tomography of cold atomic clouds could 
also be used for the observation of non-trivial density 
profiles, such as the s-wave scattering shell of two col- 
liding clouds, vortices, solitons, or the fringe pattern of 
two interfering condensates. In cases of cylindrical sym- 
metry, computerized tomography requires only a single 
absorption image. 
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